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ABSTRACT 


Multi-junction  solar  cells  are  an  emerging  technology  that  improves  the 
conversion  rate  of  solar  energy.  Indium  Gallium  Phosphide  (InGaP)  is  commonly 
used  as  the  top  cell  in  multi-junction  cells  grown  on  Germanium  (Ge)  or  Gallium 
Arsenide  (GaAs)  substrates.  To  design  more  efficient  solar  cells  using  InGaP,  it 
is  important  to  characterize  its  transport  parameters,  particularly  the  minority 
charge  carrier  mobility,  diffusion  length  and  lifetime  as  a  function  of  doping  and 
material  growth  conditions. 

In  this  work,  transport  imaging  was  performed  on  a  set  of  InGaP 
heterostructures  (with  differing  thicknesses,  doping  levels  and  minority  carrier 
types)  to  determine  their  minority  carrier  diffusion  length.  These  measurements, 
together  with  an  independent  set  of  time-resolved  photoluminescence  (TRPL) 
lifetime  data,  were  used  to  calculate  the  minority  carrier  mobility  values.  For  the 
shortest  diffusion  lengths,  experimental  limitations  were  encountered  involving 
the  finite  carrier  generation  volume.  Simulations  were  performed  to  explore  the 
potential  of  modeling  the  convolution  of  diffusion  behavior  with  a  finite  generation 
region  to  address  these  limitations. 

Transport  imaging  was  also  performed  on  a  set  of  Copper  Indium  Gallium 
Selenide  (CIGS)  materials.  Polycrystalline  CIGS  represents  an  alternative  to  the 
expensive  single-crystal  InGaP.  These  initial  experiments  identified  the 
challenges  of  applying  transport  imaging  to  polycrystalline  materials. 
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I.  INTRODUCTION 


A.  SOLAR  CELL  TECHNOLOGY 

The  expanding  use  of  solar  cell  technology  promises  a  source  of  energy 
that  is  environmentally  friendly,  renewable,  and  most  importantly,  lasting.  In  this 
era,  there  is  an  ever-increasing  power  consumption  requirement  to  fuel  the  global 
economy  and  its  operations.  While  fossil  fuel  continues  to  be  the  mainstay  of 
energy  providers,  this  energy  source  is  not  limitless  and  will  eventually  be 
depleted.  In  addition,  the  by-products  of  burning  fossil  fuel  contribute 
significantly  to  pollution,  which  subsequently  harms  the  eco-system  in  which  we 
live.  As  such,  solar  energy  can  serve  as  a  viable  option  for  replacing  fossil  fuel, 
in  part  if  not  in  whole,  with  fewer  of  its  detrimental  effects. 

While  solar  energy  presents  immense  potential,  there  are  challenges  to 
how  we  can  effectively  tap  this  energy  source.  The  first  solar  cell  was  developed 
by  Charles  Fritts  in  the  1880s.  At  its  earliest  stage,  it  could  effectively  convert 
only  less  than  1%  of  the  incident  light  [1],  As  technology  progressed,  current 
single-junction  solar  cells  are  now  able  to  attain  a  conversion  rate  of 
approximately  10  to  20%,  with  the  highest  efficiency  record  held  at  27.6%  [2], 
Conversion  efficiency  remains  an  important  issue  when  considering  the  adoption 
of  solar  cells  on  a  large  scale.  The  existing  conversion  rate  makes  solar  cells 
less  efficient  and  requires  more  panels  to  meet  a  given  requirement.  This,  in 
turn,  drives  the  implementation  cost  up,  making  it  less  appealing  as  a  substitute 
for  traditional  fossil  fuel.  The  aim  for  solar  cell  research  is  therefore  to 
simultaneously  raise  the  limit  of  conversion  efficiency  and  drive  down 
manufacturing  costs. 

Solar  cells  are  photovoltaic  devices  that  absorb  sunlight  and  convert  the 
energy  from  the  electromagnetic  spectrum  into  usable  electricity  via  the 
photovoltaic  effect.  Designed  essentially  as  semiconductors  with  p-n  junctions, 
free  charge  carriers  are  created  within  the  material  with  the  absorption  of  incident 
light.  To  achieve  high  conversion  efficiency,  these  charges  will  need  to  be 
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efficiently  transported  out  of  the  bulk  material  before  they  are  lost,  primarily  due 
to  the  generation-recombination  effect.  The  ability  of  these  charges  to  be 
transported  is  primarily  dependent  on  the  diffusion  length  of  minority  carriers. 
The  time  that  a  minority  charge  carrier  exists  before  it  recombines  is  known  as 
the  lifetime,  and  the  distance  that  the  charge  travels  in  its  lifetime,  based  on  its 
diffusivity,  is  called  the  diffusion  length.  Based  on  lifetime  and  diffusion  length, 
the  mobility  of  the  electron  or  hole  in  the  lattice  can  be  determined,  thus  allowing 
better  understanding  of  the  diffusivity  of  the  minority  charge  carriers  in  a  material. 

There  are  many  challenges  to  building  an  efficient  solar  cell.  Key 
parameters  in  characterizing  solar  cells  start  with  the  material  charge  transport 
parameters,  namely  the  charge  mobility,  lifetime  and  diffusion  length.  It  is  the 
goal  of  every  solar  cell  manufacturer  to  maintain  high  carrier  lifetime  so  as  to 
realize  long  diffusion  lengths.  Long  diffusion  lengths  ensure  that  the  charges, 
generated  from  the  incident  light,  can  effectively  travel  through  its  generation 
region  to  reach  the  depletion  region  where  they  can  be  collected.  Thus,  being 
able  to  accurately  determine  the  diffusion  length  of  any  solar  cell  material  proves 
advantageous  by  providing  insights  for  the  future  design  and  manufacture  of  a 
highly  efficient  solar  cell.  A  rapid  characterization  tool  that  can  determine  the 
diffusion  length  directly  following  material  growth  could  be  very  valuable. 

B.  PURPOSE  OF  THESIS 

The  objective  of  this  thesis  is  to  determine  the  minority  carrier  diffusion 
length  and  mobility  of  several  Indium  Gallium  Phosphide  (InGaP)  materials  for 
use  in  multi-junction  solar  cells.  This  will  be  done  through  transport  imaging,  a 
contact-free,  optical  technique  that  directly  images  the  minority  carrier  diffusion 
within  the  samples.  The  method  uses  a  scanning  electron  microscope  (SEM) 
coupled  with  an  optical  microscope  to  image  the  recombination  luminescence 
due  to  charge  motion.  From  the  images,  the  minority  carrier  diffusion  length  can 
be  easily  extracted.  Together  with  independent  carrier  lifetime  data,  determined 
by  time-resolved  photoluminescence  (TRPL),  the  carrier  mobility  can  be 
calculated. 
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The  determination  of  these  key  material  charge  parameters  will  enable  the 
full  characterization  of  the  solar  cell  materials  and  provide  a  better  understanding 
as  to  the  performance  of  the  solar  cells.  Furthermore,  direct  measurements  of 
the  behavior  of  minority  carriers  in  InGaP  are  limited  and  experimental  data  on 
the  mobility  of  these  carriers  are  also  scarce  in  existing  literature.  The  work  done 
in  this  thesis  will  provide  an  opportunity  to  investigate  and  further  determine 
experimental  values  of  such  carriers  in  the  solar  cell  samples. 

C.  MILITARY  RELEVANCE 

Power  sources  are  an  essential  part  of  any  military  operations,  especially 
in  the  current  context  where  reliance  on  machineries  and  systems  is  becoming 
more  prevalent.  In  order  to  sustain  such  operations,  long  logistical  trains  have  to 
be  established  for  the  supplying  of  sufficient  energy.  To  add  to  the  complication, 
operations  nowadays  also  tend  to  occur  far  away  from  the  home  base,  posing 
serious  challenges  to  the  current  modus  operandi  of  logistical  support  in  terms  of 
energy.  What  is  needed  is  a  new  fuel  alternative  that  is  renewable  and  can  be 
carried  by  the  front-line  troops,  making  them  self-sustainable.  Efficient  solar  cell 
devices  can  be  one  answer  to  this  problem. 

Ray  Mabus,  the  U.S.  Secretary  of  Navy,  has  commented  that  there  is 
now,  more  than  ever,  an  impetus  for  the  U.S.  military  to  look  at  renewable 
energy  [3],  This  issue  has  implications  at  both  the  strategic  and  tactical  levels. 
Strategically,  traditional  fossil  fuel  can,  most  often,  only  be  bought  from  regions 
which  are  either  volatile  or  have  the  potential  to  become  so.  Any  sovereign 
state’s  reliance  on  such  fuel  sources  creates  an  extreme  vulnerability  to  its 
overall  national  security.  Tactically,  it  is  highly  costly  to  supply  such  fossil  fuel  to 
front-line  troops.  The  cost  is  not  only  in  a  monetary  sense,  but  also  involves  the 
lives  of  soldiers  or  Marines  lost  in  trying  to  bring  these  fuels  into  the  theater  of 
operation.  By  adopting  solar  energy,  a  state  is  free  of  geo-political  constraints  at 
the  strategic  level.  At  the  tactical  level,  ground  commanders  can  be  assured  of 
self-sustainability  without  having  to  be  duly  concerned  with  or  crippled  by  the  lack 
of  logistical  support  for  energy. 
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For  Singapore,  the  energy  issue  is  also  of  major  concern  to  both  its 
political  and  military  leaders.  Being  a  small  island  state  without  any  natural 
resources,  it  can  only  look  to  foreign  parties  for  the  acquisition  for  such  energy 
sources.  As  discussed  earlier,  this  results  in  a  national-level  vulnerability  that 
has  to  be  addressed.  While  Singapore  enjoys  good  ties  with  countries  both 
regionally  and  globally,  being  self-sufficient  remains  a  priority  in  the  city  country’s 
development  plan.  Solar  energy  will  prove  a  hopeful  alternative  to  meet  that  end. 

Another  area  in  which  solar  energy  proves  useful  is  the  space  arena. 
More  and  more  states  are  embarking  on  space  programs,  which  seek  to  put 
satellites  up  in  space.  For  the  military,  such  satellites  play  in  an  important  role  for 
communications  as  well  as  intelligence,  surveillance  and  reconnaissance  (ISR) 
purposes.  These  satellites  are  usually  in  geo-synchronous  orbits  and  are 
constantly  exposed  to  the  sun.  Multi-junction  solar  cells  provide  a  feasible 
means  of  power  generation  that  is  both  stable  and  cost  effective. 

D.  THESIS  OVERVIEW 

Chapter  I  begins  by  giving  a  brief  introduction  to  solar  cell  technology  and 
its  current  state  of  development.  It  also  outlines  the  purpose  of  this  thesis  and 
the  work  to  be  done.  Finally,  it  evaluates  the  relevance  of  solar  cell  technology 
as  applied  to  the  military. 

Chapter  II  presents  background  information  on  multi-junction  solar  cells 
using  InGaP  and  how  they  work.  The  basic  concepts  of  semiconductor 
parameters,  such  as  minority  carrier  mobility,  lifetime  and  diffusion  length,  are 
reviewed  in  greater  depth. 

Chapter  III  explains  the  technique  of  transport  imaging.  It  also  describes 
the  setup  of  the  laboratory  equipment  as  well  as  the  experimental  procedures  for 
applying  the  technique  to  the  samples. 

Chapter  IV  summarizes  the  results  obtained  from  the  experiments  and 
provides  the  analysis  of  the  results.  It  also  highlights  the  experimental  limitations 
encountered  involving  finite  generation  volume  for  charge  carriers  for  the 
samples  with  short  diffusion  lengths. 
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Chapter  V  explores  a  deconvolution  method  to  overcome  the  limitations 
encountered  for  finite  charge  carrier  generation  volume.  The  deconvolved 
results  are  evaluated  and  compared  against  the  initial  results. 

Chapter  VI  introduces  the  Copper  Indium  Gallium  Selenide  (CIGS) 
material  system  as  an  alternate  type  of  solar  cell.  The  results  of  an  initial 
investigation  of  the  transport  imaging  measurement  of  CIGS  are  presented. 
Some  of  the  challenges  of  applying  the  transport  imaging  technique  on  CIGS  are 
also  identified. 

Finally,  Chapter  VII  concludes  the  report  with  a  summary  of  the  work  done 
and  gives  the  direction  for  which  future  research  relating  to  the  thesis  can  be 
carried  out. 
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II.  BACKGROUND 


A.  MAJORITY  AND  MINORITY  CHARGE  CARRIERS  IN 

SEMICONDUCTORS 

In  all  crystalline  structures,  the  interactions  of  atoms  or  molecules  within 
the  crystal  lattice  result  in  discrete  energy  levels,  which  combine  to  form  energy 
bands.  In  semiconductors,  the  lower  occupied  energy  bands  are  known  as  the 
valence  bands  while  the  higher  ones  are  known  as  the  conduction  bands.  In 
between  the  conduction  and  valence  bands  is  a  range  of  forbidden  energy  levels 
called  the  band-gap.  Within  the  valence  band,  the  electrons  are  bound  by 
energy  bonds  to  their  parent  atoms  in  the  crystal  lattice. 

When  an  electron  gains  enough  energy,  either  through  thermal  excitation 
or  photon  absorption,  it  is  possible  for  the  electron  to  break  free  of  its  chemical 
bonds  and  make  an  upward  transition  to  the  conduction  band.  The  condition  for 
this  to  happen  is  that  the  energy  gained,  A E ,  must  be  greater  than  or  equal  to  the 
band-gap  energy,  e  ,  i.e.,  the  energy  needed  for  it  to  cross  the  forbidden  band- 

gap  (  ae  >  Eg  ).  The  electron  is  then  said  to  have  make  a  “quantized  jump.”  After 

the  transition,  there  is  a  vacancy  left  behind  in  the  valence  band  by  the  electron, 
called  a  hole.  This  is  the  generation  of  an  electron-hole  pair  (Figure  1 ).  Both  the 
electron,  negatively  charged,  and  the  hole,  positively  charged,  can  then  move 
freely  within  their  respective  bands,  especially  under  the  influence  of  an  applied 
electric  field.  The  term  “charge  carrier”  is  used  to  refer  to  both  the  electrons  and 
holes. 
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Figure  1 .  Generation  of  Electron-Hole  Pair  (After  [4]) 

An  intrinsic,  or  pure,  semiconductor  will  have  an  equal  number  of 
electrons  and  holes  across  the  band-gap  under  steady  state  thermal  equilibrium. 
The  reason  for  this  is  that  electrons  and  holes  are  always  created  in  pairs.  The 
charge  carrier  (electron  or  hole)  concentration  is  governed  by  the  following 
equation: 

noPo=n t  (Eq.  1) 

where  n0  is  the  electron  charge  carrier  concentration,  p0  is  the  hole  charge 
carrier  concentration  and  n ;.  is  the  average  intrinsic  charge  carrier  concentration. 

To  improve  and  control  the  conductivity  of  the  semiconductor,  the  doping 
process,  which  is  the  intentional  introduction  of  impurities,  is  often  performed  on 
the  material.  In  this  case,  the  semiconductor  now  becomes  extrinsic.  Doping 
can  either  be  done  with  donor  or  acceptor  dopants. 

Donor  dopants  (e.g.,  Group  V  elements  for  doping  of  Si)  have  one  extra 
electron  than  is  required  for  bonding  with  the  surrounding  atoms  in  the  crystal 
lattice.  As  such,  an  electron  is  deemed  have  been  ‘donated’  and  is  free  to  move 
in  the  conduction  band.  On  a  macro  scale,  the  overall  concentration  of  electrons 
will  now  be  greater  than  the  concentration  of  holes  ( n0>p0 )  at  room 

temperature.  The  electrons,  being  of  higher  concentration,  are  called  the 
majority  charge  carriers,  and  the  holes,  being  of  lower  concentration,  are  called 
the  minority  charge  carriers.  The  semiconductor  material  is  now  referred  to  as 
an  n-type  material  (Figure  2). 
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Figure  2.  Example  of  a  N-Type  Material  (Silicon  Doped  with  Antimony)  (From 

[5]) 

Acceptor  dopants  (e.g.,  Group  III  elements  for  doping  of  Si)  have  one  less 
electron  than  is  needed  for  bonding.  They  will  then  readily  ‘accept’  a  nearby 
electron,  creating  an  extra  hole  in  the  process  (Figure  3).  Again,  on  a  macro 
level,  this  causes  the  concentration  of  holes  to  be  greater  than  the  concentration 
of  electrons  ( p0  >n0)  at  room  temperature.  Therefore,  the  holes  are  now  the 

majority  charge  carriers  while  the  electrons  are  the  minority  charge  carriers,  and 
the  semiconductor  material  is  known  as  a  p-type  material. 
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Figure  3.  Example  of  a  P-Type  Material  (Silicon  Doped  with  Boron)  (From  [5]) 


The  behavior  of  minority  charge  carriers,  rather  than  the  majority  charge 
carriers,  is  more  often  the  primary  concern  of  semiconductor  manufacturers  or 
solar  cell  designers.  When  exposed  to  photons  with  energy^^,  both  types  of 

charge  carriers  are  generated  (electron-hole  pair  generation).  For  both  n-type 

and  p-type  materials,  the  change  in  concentration  for  the  majority  carrier  due  to 
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the  excitation  is  generally  very  small  because  of  its  large  existing  carrier 
population.  On  the  other  hand,  the  fractional  change  in  concentration  for  the 
minority  carrier  is  considerably  more  significant  in  comparison.  As  a  result, 
minority  charge  carriers  play  a  more  critical  role  in  determining  the  non¬ 
equilibrium  electrical  properties  of  the  material,  making  it  an  important  parameter 
for  consideration. 

B.  CHARGE  CARRIER  DIFFUSION  LENGTH 

Diffusion,  in  general,  is  the  migration  of  particles  from  a  region  of  higher 
concentration  to  a  region  of  lower  concentration.  It  is  akin  to  the  spread  of  a  drop 
of  ink  in  a  glass  of  water  or  the  molecular  spreading  of  a  gas  in  a  container  [6], 
When  a  semiconductor  is  exposed  to  illumination,  generation  of  electron-hole 
pairs  occurs.  This  is  usually  localized,  resulting  in  a  region  which  suddenly 
experiences  a  higher  concentration  of  charge  carriers  than  the  surrounding 
material.  This  difference  in  concentration  creates  a  concentration  gradient  and 
causes  the  charge  carriers  to  diffuse  into  the  immediate  vicinity.  Figure  4  shows 
the  migration  of  the  charge  carriers  away  from  the  center  point  of  illumination, 
smoothing  out  the  concentration  gradient  over  time.  For  majority  carriers,  there 
only  need  be  a  small  change  to  the  local  concentration  to  accommodate  the  flux 
of  the  diffusing  carriers,  whereas  there  is  a  more  significant  change  for  the 
minority  carriers.  Once  again,  the  importance  of  understanding  minority  carrier 
diffusion  in  semiconductors  is  highlighted. 
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Figure  4.  Charge  Carrier  Diffusion.  The  top  Portion  Shows  the  Generation  of 
Electron-Hole  Pairs  due  to  Incident  Light.  The  Bottom  Portion  Shows  the 
Broadening  of  the  Carrier  Distribution  Over  Time  (From  [6]) 


The  diffusivity  of  charge  carriers  within  a  material  is  described  by  the 
diffusion  coefficient,  D .  This  coefficient  is  governed  by  the  Einstein  relation  as 
given  below: 


D 


/AT 


(Eg.  2) 


e 

where  /u  is  the  charge  carrier  mobility,  k  is  Boltzmann’s  constant,  and  T  is  the 
temperature.  An  important  thing  to  note  is  that  the  mobility  value  for  electrons 
and  holes  are  different  even  if  within  the  same  material.  This  is  due  to  their 
difference  in  size  and  mass  and  also  as  a  result  of  their  transport  in  different 
bands.  Therefore,  the  diffusion  coefficients  for  electrons  and  holes  are  also 
different. 

The  diffusion  length  is  defined  as  the  average  length  a  charge  carrier,  be  it 
electron  or  hole,  is  able  to  diffuse  before  it  recombines  and  returns  to  its 
equilibrium  state.  The  diffusion  length,  z  can  be  calculated  by  the  following 


equation: 
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(Eq.  3) 


where  is  D  the  diffusion  coefficient  as  defined  earlier,  and  r  is  generally  the 
minority  charge  carrier  lifetime. 

C.  CHARGE  CARRIER  MOBILITY 

The  mobility  is  a  measure  of  the  relative  ease  with  which  a  charge  carrier 
can  move  within  a  material.  In  general,  the  mobility,  // ,  can  be  found  by: 

PT 

M  =  ^  (Eq.  4) 

m 

where  e  is  the  electron  charge,  zs  is  the  scattering  lifetime  and  m*  is  the 

effective  mass.  The  difference  in  mobility  between  electrons  and  holes  can  be 
explained  by  the  fact  that  electrons  and  holes  have  different  effective  masses, 
resulting  in  different  mobility  values. 

Mobility  of  minority  charge  carriers  is  traditionally  determined  by  the 
Haynes-Shockley  technique  [7],  In  the  experiment,  a  pulse  of  charge  carriers  is 
created  at  one  end  of  the  semiconductors  (usually  done  by  having  light  incident 
on  it).  Under  an  applied  electric  field,  the  charge  carriers  will  be  made  to  move 
across  the  length  of  the  semiconductor  to  the  other  end.  Using  an  oscilloscope, 
the  transit  time  can  be  found.  The  velocity  of  the  carrier  can  then  be  easily 
calculated  since  the  distance  traveled  (which  is  equal  to  the  length  of  the 
semiconductor)  and  the  transit  time  is  known.  The  mobility  is  defined  as  the 
proportionality  constant  between  the  applied  electric  field,  E,  and  the  drift 
velocity,  v: 

(Eq.  5) 

E 

The  disadvantages  to  this  technique  are  that  contacts  have  to  be  added  to 
the  semiconductor  to  induce  the  needed  electric  field  and  that  it  averages  over 
the  full  length  of  the  sample. 
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The  mobility  can  be  influenced  by  both  temperature  and  the  doping  level 
of  the  semiconductor.  As  the  charge  carriers  move  through  the  material,  they 
can  be  scattered  by  the  atoms  vibrating  about  their  fixed  positions  in  the  crystal 
lattice.  When  temperature  increases,  the  atoms  gain  more  kinetic  energy  and 
the  vibrational  displacement  increases.  This  causes  a  higher  possibility  of 
collision  between  the  charge  carriers  and  the  vibrating  atoms,  affecting  the 
mobility.  With  regard  to  doping,  charged  ions  are  introduced  to  the  lattice  in 
place  of  neutral  atoms.  The  charge  of  the  ions  has  either  an  attractive  or 
repulsive  effect  on  the  moving  charge  carriers,  acting  to  deflect  and  scatter  them, 
thereby  affecting  the  mobility.  Hence,  the  mobility  is,  in  certain  temperature 
regimes,  inversely  proportional  to  the  doping  level. 

D.  CHARGE  CARRIER  LIFETIME 

The  carrier  lifetime  is  the  mean  time  that  a  carrier  exists  before  it 
recombines  and  disappears.  A  non-equilibrium  charge  carrier  does  not  exist 
forever.  As  it  travels  through  the  bulk  material,  it  will  eventually  encounter  an 
opposite  charge  carrier  (e.g.,  an  electron  meeting  a  hole)  and  recombine  with  it. 
Therefore,  for  the  minority  carriers,  the  lifetime  is  commonly  inversely 
proportional  to  the  population  of  the  majority  carriers: 


where  B  is  the  recombination  coefficient  and  N  is  the  majority  carrier 
population. 

Charge  carrier  lifetime  can  be  determined  through  the  time-resolved 
photoluminescence  (TRPL)  technique  [8],  A  variety  of  configurations  for  TRPL 
exist,  but  generally,  the  semiconductor  sample  is  excited  with  a  pulsed  external 
illumination  source  and  the  subsequent  photoluminescence  (PL)  against  time 
profile  of  the  semiconductor  is  observed  and  recorded.  The  lifetime  can  be 
evaluated  from  the  decay  in  PL  intensity  by: 

(Eq.7) 
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where  /0is  the  initial  intensity  while  /  is  the  instantaneous  intensity  at  time  t . 


E.  LUMINESCENCE  IN  SEMICONDUCTORS 

Luminescence  in  semiconductors  arises  when  the  recombination  process 
between  an  electron  and  a  hole  results  in  the  emission  of  a  photon,  releasing 
energy  concurrently  (Figure  5).  As  explained  earlier,  an  excited  electron,  usually 
created  through  the  absorption  of  incident  photon  energy,  resides  in  the 
conduction  band.  Eventually,  it  will  meet  with  a  hole  and  recombine  with  it, 
making  a  transition  from  its  higher  energy  state  into  the  lower,  but  more  stable, 
valence  state.  The  recombination  affords  the  electron  a  means  to  liberate  its 
extra  energy  in  the  form  of  a  photon  emission.  Hence,  luminescence  can  be 
observed.  The  energy  of  the  emitted  photon  is  approximately  equal  to  the  band- 
gap  energy  of  the  semiconductor  and  the  photon  wavelength  can  be  found  using 
the  following  equation: 

1.24 

E*=-jT  (Eq-8) 

where  e  is  the  band-gap  energy  (in  eV),  and  the  A  is  the  wavelength  (in  pm). 


Electron  Energy,  E 


Figure  5.  Photon  Absorption  and  Recombination  (From  [9]) 


Each  semiconductor  material  is  uniquely  characterized  by  its  band-gap 
energy.  Therefore,  it  is  possible  to  determine  the  material  of  a  semiconductor  by 
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using  spectroscopy  to  analyze  the  wavelength  of  the  emitted  light.  The  band-gap 
energy  for  semiconductors  is  typically  in  the  range  of  0.1  eV  to  5  eV. 


F.  MULTI-JUNCTION  SOLAR  CELLS 

A  solar  cell  works  by  converting  energy  from  the  sun  into  usable 
electricity.  When  the  solar  cell  is  exposed  to  sunlight,  the  semiconductor  material 
in  the  device  absorbs  energy  from  the  solar  radiation.  This  energy  is  used  to 
create  charge  carriers,  which  are  free  to  flow  within  the  solar  cell.  Due  to  the  p-n 
junctions,  there  is  a  built-in  electric  field,  which  causes  the  charge  carriers  to 
move  in  a  particular  direction,  thereby  creating  an  electric  current.  By  placing 
metal  contacts  at  the  top  and  bottom  of  the  solar  cell,  the  current  is  allowed  to 
flow  out  of  the  device  for  use  externally  [1 0], 

Energy  from  the  sun  consists  of  photons  with  a  multitude  of  wavelengths 
corresponding  to  the  electromagnetic  solar  spectrum.  A  solar  cell,  no  matter  how 
efficient,  can  only  work  to  convert  a  portion  of  the  sun’s  energy.  The  sole  factor 
governing  the  energy  that  can  be  absorbed  is  the  material’s  band-gap 
characteristic.  In  the  case  of  a  single-junction  solar  cell,  there  is  only  one  band- 
gap  threshold.  Photons  of  energy  lower  than  the  band-gap  threshold  pass 
through  the  solar  cell  without  being  absorbed,  since  they  have  insufficient  energy 
to  excite  the  electrons  into  the  conduction  band.  For  photons  with  energy  higher 
than  the  threshold,  a  portion  of  the  energy  is  used  for  the  electron  to  overcome 
the  band-gap  and  cross  to  the  conduction  band,  while  any  excess  energy  is 
converted  to  heat  and  is  dissipated  throughout  the  bulk  material.  This  is 
unfavorable  as  it  increases  the  lattice  vibration  in  the  semiconductor,  leading  to 
shorter  lifetime  of  the  charge  carriers  and  therefore,  shorter  diffusion  length.  This 
loss  of  energy  also  limits  the  overall  efficiency  of  the  device. 

To  absorb  more  photons  of  different  energies  and  improve  the  conversion 
efficiency,  multi-junction  solar  cells  are  designed.  This  is  achieved  by  stacking 
various  types  of  solar  cells  into  a  single  device.  Each  cell  has  a  distinct  band- 
gap  energy  and  is  thus  tuned  to  absorb  photons  with  energy  greater  than  a 

specific  cut-off  wavelength.  In  general,  the  solar  cells  are  stacked  in  decreasing 
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order  of  band-gap  threshold.  The  top  cell  will  be  able  to  absorb  photons  with  the 
highest  energies,  while  the  lower  energy  photons  will  pass  through  and  be 
subsequently  absorbed  by  the  lower  cells,  each  with  a  lower  band-gap  threshold 
than  the  one  before  (Figure  6).  Together,  the  multi-junction  cell  can  now  absorb 
more  of  the  solar  radiation  while,  at  the  same  time,  minimizing  energy  which  is 
lost  as  heat  [11],  Ideally,  a  multi-junction  solar  will  be  made  up  of  hundreds  of 
layers,  each  able  to  absorb  a  certain  part  of  the  electromagnetic  spectrum. 
However,  due  to  technological  limitations,  current  multi-junction  solar  cells  are 
designed  with  no  more  than  just  a  few  layers  (triple-junction  solar  cell  is  the  most 
common  type  at  the  time  of  writing)  [1 2], 


Figure  6.  Schematic  of  a  Triple-Junction  Solar  Cell.  The  Top  Cell  Absorbs 
Photons  with  the  Highest  Energy,  £  while  the  Subsequent  Cells  Absorb 

Photons  with  the  Lower  Energies,^,,  and  Eg3  (From  [12]) 

lno.5Gao.5P  is  a  III  -  V  ternary  semiconductor  made  up  of  the  elements 
indium,  gallium  and  phosphorous.  It  is  often  used  as  the  top  layer  in  the 
fabrication  of  multi-junction  solar  cells  due  to  its  large  band-gap  characteristic 
and  its  ability  to  be  grown  lattice-matched  to  Gallium  Arsenide  (GaAs).  It  has  a 
near-perfect  lattice  match  with  GaAs  (lattice  mismatch  is  generally  at  about  0.2% 
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[13]),  another  preferred  semiconductor  for  solar  cell  manufacturing.  This  enables 
InGaP  to  be  easily  grown  on  GaAs,  and  the  two  are  commonly  used  as  tandem 
solar  cells. 
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III.  TRANSPORT  IMAGING 


A.  DIRECT  IMAGING  OF  FREE  CARRIER  DIFFUSION 

Transport  imaging  is  an  optical  imaging  technique  that  seeks  to  directly 
measure  the  diffusion  length  of  minority  carriers  by  imaging  the  spatial 
dependence  of  recombination  luminescence  of  a  material.  In  this  technique,  an 
external  source  will  be  used  to  locally  excite  the  material,  creating  electron-hole 
pairs.  This  generates  a  carrier  gradient  at  the  point  of  excitation  and  induces  the 
subsequent  diffusion  of  the  minority  charge  carriers.  The  image  of  the 
luminescence  distribution,  as  a  result  of  the  radiative  recombination  of  the 
carriers  as  they  diffuse,  is  captured  using  a  charge-coupled  device  (CCD) 
camera.  Unlike  the  traditional  cathodoluminescence  (CL)  method,  the  spatial 
information  of  the  luminescence  distribution  is  retained  [14],  allowing  for  the 
extraction  of  the  diffusion  length  and  the  subsequent  calculation  of  the  minority 
charge  carrier  mobility. 

One  advantage  of  transport  imaging  over  the  conventional  Haynes- 
Shockley  technique  is  that  contacts  do  not  need  to  be  added  to  the  sample  in 
order  to  induce  an  electric  field  necessary  to  move  the  minority  charge  carriers. 
This  is  also  the  advantage  over  electron-beam-induced  current  (EBIC) 
techniques,  which  require  a  rectifying  contact  to  collect  electron-beam  induced 
charge.  In  addition,  it  can  be  performed  with  much  higher  spatial  resolution. 

B.  LABORATORY  EQUIPMENT 

1.  SEM 

In  this  experiment,  the  external  excitation  source  will  consist  of  an  electron 
beam  generated  in  a  SEM.  The  SEM  used  is  the  JEOL  840A  model  and  it  is 
fitted  with  a  modified  stage  from  Oxford  Instruments  to  allow  the  intake  of  helium 
gas  for  the  cooling  of  the  sample  stage  (Figure  7).  This  modification  enables 
measurement  of  the  sample  to  be  conducted  at  a  range  of  temperatures — from 
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room  temperature  of  300K  to  a  low  temperature  of  5K.  There  is  also  an  optical 
microscope  attached  to  collect  the  luminescence  from  the  samples. 


Figure  7.  JEOL  840A  SEM  with  Modified  Variable  Temperature  Stage 

The  SEM  can  be  adjusted  to  operate  in  three  separate  modes — spot 
mode,  line  mode  and  picture  mode  (Figure  8).  In  spot  mode,  the  electron  beam 
is  held  fixed  at  a  specified  spot,  allowing  charge  carriers  to  be  generated  from 
that  quasi-point  source.  In  line  mode,  the  electron  beam  is  scanned  along  a 
“line”  and  charges  are  generated  along  the  beam  path.  In  picture  mode,  the 
beam  scans  over  a  specified  region,  allowing  an  area  “picture”  of  the  material 
imaged. 


Figure  8.  SEM  Imaging  Modes:  (a)  Spot  Mode,  (b)  Line  Mode  and  (c)  Picture 

Mode 
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2.  Optical  Detector 

The  optical  detector  connected  to  the  optical  microscope  for  luminescence 
imaging  is  an  Apogee  CCD  camera  (Figure  9).  It  has  a  2184  x  1472  pixel  array. 
When  coupled  to  our  optical  microscope,  the  resolution  is  0.4  pm  per  pixel.  In  an 
unfiltered  setting,  the  camera  can  efficiently  collect  data  for  light  of  wavelengths 
in  the  range  of  400  to  1100  nm.  To  minimize  the  effect  of  thermal  noise,  the 
camera  can  also  be  cooled  to  a  temperature  of  -20  °C  prior  to  operation. 


Figure  9.  Apogee  CCD  Camera 


Figure  10  shows  the  overall  operation  of  transport  imaging  measurement 
system  as  well  as  a  schematic  of  the  components  used. 


Electron  Beam 


Optical  Microscope 


Image 

Acquisition  and 
Display 


Figure  10.  Transport  Imaging  Operation  and  Components  (After  [15]) 
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C.  DATA  COLLECTION 

For  transport  imaging,  both  the  spot  mode  and  line  mode  can  be  used. 
However,  the  line  mode  was  predominantly  used  for  these  experiments  due  to 
the  relative  ease  of  the  extraction  of  diffusion  length  data  as  compared  to  that 
when  spot  mode  is  used.  The  reason  for  this  is  further  discussed  in  Section  D. 

For  each  sample,  two  separate  images — the  “luminescence”  image  and 
the  “dark”  image — were  taken  (Figure  1 1).  The  first  is  the  luminescence  image, 
which  shows  the  steady  state  luminescence  distribution  of  the  sample  when 
exposed  to  the  SEM’s  electron  beam.  The  second  is  the  dark  image  of  the 
sample  when  the  electron  beam  has  been  turned  off  or  blocked.  This  image, 
though  appearing  totally  black  to  the  naked  eye,  actually  captures  the 
background  CCD  response  due  to  thermal  generation  in  the  detector  and  any 
background  light  in  the  chamber.  During  data  processing,  the  luminescence 
count  from  the  dark  image  is  subtracted  from  the  count  of  the  luminescence 
image  to  acquire  a  more  accurate  luminescence  distribution,  which  is  due  solely 
to  the  material  itself.  To  obtain  a  higher  level  of  accuracy,  both  images  are  taken 
over  the  same  region  of  the  sample  and  exposed  for  the  same  duration  of  time. 


Figure  1 1 .  Data  Collection  for  a  10s  Exposure:  (a)  the  Luminescence  Image  and 
(b)  the  Dark  Image  [712  pixels  (285  pm)  x  712  pixels  (285  pm)] 
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Another  important  factor  to  ensure  the  proper  capturing  of  the 
luminescence  distribution  is  to  have  a  good  signal-to-noise  ratio  (SNR).  For  the 
experiment,  the  camera  is  exposed  to  the  sample  for  a  period  of  time  (which  will 
vary  between  samples  due  to  their  different  luminescence  intensity)  such  that  the 
peak  intensity  reaches  a  count  of  approximately  25,000  or  higher.  With  the 
background  count  at  about  200,  a  good  SNR  of  close  to  50dB  can  be  achieved. 
Figure  12  shows  the  horizontal  intensity  profile  of  a  sample  luminescence  image 
in  line  mode. 


Figure  12.  Intensity  Profile  of  a  Luminescence  Image  (Extracted  Perpendicular 

to  the  Excitation  Line) 

D.  DIFFUSION  LENGTH  EXTRACTION 
1.  2D  Diffusion  Case 

The  extraction  method  for  the  diffusion  length  depends  on  the  imaging 
mode  with  which  the  sample  was  taken.  In  the  case  of  a  spot  mode,  the  minority 
charge  carriers  are  allowed  to  diffuse  from  a  fixed  generation  point.  Assuming 
that  the  sample  is  a  thin  layer  and  that  the  diffusion  is  thus  confined  to  a  plane, 
the  charge  carriers  will  follow  a  two-dimensional  (2D)  pattern  of  diffusion  (Figure 
13). 
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Electron  Beam 


Side  View 


Top-down  View 


Figure  13.  2D  Diffusion  of  Charge  Carriers 

This  diffusion  can  be  modeled  as  a  zeroth-order  modified  Bessel  function 
of  the  second  kind  [14]  and  the  continuity  equation  involving  the  minority  carrier 
distribution  and  related  recombination  luminescence,  /  ,  can  be  expressed  as: 

0  =  DW2I  -  —I  +  S{r)  (Eq.9) 

T 

where  D  is  the  minority  charge  carrier  diffusivity,  r  is  the  minority  charge  carrier 
lifetime  and  is  S(r )  the  source  function. 

The  luminescence,  / ,  resulting  from  a  point  source,  S(r),  can  then  be 
defined  by  solving  the  differential  equation  given  in  Equation  9  and  subsequently 
approximated  by: 

W  =  (Eq.  10) 

2  7tD 

where  g  is  the  amplitude  of  the  source,  K0  is  the  zeroth-order  modified  Bessel 
function  of  the  second  kind,  r  is  the  radius  of  diffusion  and  Ldiff  is  the  diffusion 
length.  The  diffusion  length,  Ldiff,  is  determined  from  a  least  squares  analysis  for 
variations  in  Ldiff . 
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2. 


ID  Diffusion  Case 


In  the  case  of  a  line  mode,  the  electron  beam  scans  up  and  down  the 
sample  along  a  straight  path,  causing  the  charge  carriers  to  diffuse  from  a  line 
source.  Since  only  the  exitance  condition  is  different  for  the  charge  carriers 
diffusing  at  either  ends  of  the  line,  the  charge  carriers  in  the  middle  (as  depicted 
by  the  yellow  box  in  Figure  14)  can  be  considered  to  be  following  a  one¬ 
dimensional  (1 D)  diffusion  model. 


Figure  14.  ID  Diffusion  of  Charge  Carriers 

The  1 D  minority  charge  distribution  falls  off  exponentially  as  a  function  of 
the  distance  and  the  luminescence,  /  ,  can  be  expressed  as: 

I  =  I0eLdV  (Eq.  11) 

where  I0  is  the  peak  intensity,  x  is  the  distance  of  diffusion  from  the  line  and 
Ldiff  is  the  diffusion  length. 
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After  normalizing  luminescence,  the  diffusion  length  can  be  calculated  by: 

Ax 

(Eg.  12) 


off 


A  In  A 


where  lN  is  the  normalized  luminescence  (i.e.  IA 


/~4,n 

^0  —  I  mm 


)■ 


If  the  semi-logarithmic  graph  of  the  normalized  luminescence  intensity,  I N , 

against  position,  x,  is  plotted,  then  the  diffusion  length,  Ldiff ,  can  be  directly 

determined  as  the  inverse  slope  value  of  the  graph  (Figure  15).  This  is  known  as 
the  1/slope  technique. 
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(a)  Plot  of  Normalized  Intensity  against 
Position 


(b)  Semi-Logarithmic  Plot  of 
Normalized  Intensity  against  Position 


Position  (microns)  Position  (microns) 

(c)  Linear  Regression  Lines  for  (Left  and  Right)  Slope  Extraction 


Figure  15.  1/Slope  Technique  (From  [15]) 
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IV.  EXPERIMENTAL  RESULTS  OF  MEASURING  MINORITY 
CARRIER  DIFFUSION  LENGTH  IN  INDIUM  GALLIUM  PHOSPHIDE 


A.  SAMPLES 

For  this  experiment,  four  sets  of  InGaP  heterostructures  were  provided 
through  collaboration  with  other  researchers  in  the  solar  cell  field.  Three  of  the 
sets  had  p-type  InGaP  material  with  two  different  doping  levels  while  the  last  set 
had  n-type  InGaP  material.  Within  each  set  of  the  heterostructures,  there  were 
three  samples  with  varying  thicknesses  for  the  InGaP  layer  at  0.3  pm,  0.5  pm, 
and  1.0  pm,  respectively.  Therefore,  there  were  a  total  of  twelve  heterostructure 
cells  for  diffusion  length  measurement. 

1 .  Set  01 :  p-type  InGaP  (Higher  Doping) 

Set  01  samples  consisted  of  Alo.25Gao.25lrio.5P  barrier  layers  with  a  middle 
p-type  lno.5Gao.5P  layer  with  a  doping  level  of  1  x  1017  cm'3  (Figure  16).  The 
three  samples  within  this  set  had  different  thicknesses  of  0.3  pm,  0.5  pm,  and  1 .0 
pm  for  the  InGaP  layer.  These  samples  will  be  labeled  as  sample  1A,  sample 
IB,  and  sample  1C,  respectively. 


p-Aln  7qGan  7qlnn  qnP  (1  x  1017  cm  3,  50nm) 

P'Gao  so*no  5op  (1  x  !017  cm'3'  varying  thickness) 

p-AI0  25Ga0  25,no.5op  U  x  1017  cm  3>  50nm) 
GaAs  Substrate  (undoped) 


Figure  16.  Set  01  Samples  with  1  x  1 017  cm"3  P-Type  InGaP  Layer 

2.  Set  02:  p-type  InGaP  (Lower  Doping) 

Set  02  samples  consisted  of  Alo.25Gao.25lno.5P  barrier  layers  with  a  middle 
p-type  lno.5Gao.5P  layer  with  a  doping  level  of  5  x  1016  cm"3  (Figure  17).  The 
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three  samples  within  this  set  had  different  thicknesses  of  0.3  pm,  0.5  pm,  and  1 .0 
pm  for  the  InGaP  layer.  These  samples  will  be  labeled  as  sample  2A,  sample 
2B,  and  sample  2C,  respectively. 


p-AI0  25Ga0  25lno.5oP  (5  x  1016  cm  3,  50nm) 

p-Ga0  50ln0  50P  (5  x  1016  cm  3,  varying  thickness) 

p-Alp  25^ao  25^0.50^  (5  x  1016  cnrr3,  200nm) 
GaAs  Substrate  (undoped) 


Figure  17.  Set  02  Samples  with  5  x  1016  cm'3  P-Type  InGaP  Layer 

3.  Set  03:  p-type  InGaP  (Lower  Doping) 

Set  03  samples  consisted  of  Alo.25Gao.25lrio.5P  barrier  layers  interfaced  with 
a  middle  p-type  lno.5Gao.5P  layer  with  a  doping  level  of  5  x  1016  cm'3  (Figure  18). 
This  doping  level  was  the  same  as  that  of  those  in  Set  02.  The  three  samples 
within  this  set  had  different  thicknesses  of  0.3  pm,  0.5  pm,  and  1.0  pm  for  the 
InGaP  layer.  These  samples  will  be  labeled  as  sample  3A,  sample  3B,  and 
sample  3C,  respectively. 


P'Alo.25Gao.25,no.5op  (5  x  1016  cm  3,  50nm) 

P"Gao.so,no.5op  (5  x  1016  cm  3,  varying  thickness) 

P Alo.25Gao.25,no.5op  (5  x  1016  cm'3,  200nm) 
GaAs  Substrate  (undoped) 


Figure  18.  Set  03  Samples  with  5  x  1016  cm'3  P-Type  InGaP  Layer 
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4.  Set  04:  n-type  InGaP 

Set  04  samples  were  the  only  n-type  lno.5Gao.5P  heterostructures,  with  the 
highest  doping  level  of  8  x  1017  cm'3  (Figure  19).  The  three  samples  within  this 
set  had  different  thicknesses  of  0.3  pm,  0.5pm,  and  1.0  pm  each,  and  will  be 
labeled  as  sample  4A,  sample  4B,  and  sample  4C,  respectively. 


n'Alo.5o'no.5oP  (1  x  1017  cm‘3/  50nm) 
n-Ga0  50ln0  50P  (8  x  1017  cm  3,  varying  thickness) 
n-Al0.5o* no 50^  U  x  1017  cnr3,  200nm) 

GaAs  Substrate  (undoped) 

Figure  19.  Set  04  Samples  with  8  x  1017  cm'3  N-Type  InGaP  Layer 
A  summary  of  the  twelve  samples  is  given  in  Table  1 . 


Set# 

Type 

Sample  Name 

Thickness 

(pm) 

1 

p-InGaP 

(Higher  doping  -  1  x  1017  cm'3) 

1 A 

0.3 

IB 

0.5 

1C 

1.0 

2 

p-InGaP 

(Lower  doping  -  5  x  1 016  cm'3) 

2A 

0.3 

2B 

0.5 

2C 

1.0 

3 

p-InGaP 

(Lower  doping  -  5  x  1016  cm'3) 

3A 

0.3 

3B 

0.5 

3C 

1.0 

4 

n-InGaP 

(Doping  -  8  x  1017  cm'3) 

4A 

0.3 

4B 

0.5 

4C 

1.0 

Table  1 .  InGaP  Heterostructure  Sample  Sets 
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B.  CATHODOLUMINESCENCE  SPECTRA  OF  SAMPLES 


Prior  to  the  transport  imaging  experiment,  CL  spectroscopy  was 
performed  to  determine  the  radiative  spectrum  of  the  samples.  Figure  20  shows 
the  room  temperature  CL  spectra  results  of  the  1  pm  thick  samples  from  each 
set.  From  the  graph,  it  is  observed  that  photoluminescence  occurred  at  only  one 
peak  wavelength  for  each  of  the  sample  sets.  These  luminescence  peaked  at  a 
wavelength  of  approximately  650  nm,  which  corresponds  to  the  InGaP  bandgap 
of  about  1.9  eV  [16].  The  slight  shift  in  wavelength  is  likely  a  result  of  the 
different  doping  levels  of  the  sample  sets.  This  ensures  that  the 
photoluminescence  to  be  used  for  the  transport  imaging  originated  purely  from 
the  InGaP  layer  and  that  there  was  no  higher  energy  interference  from  the  barrier 
layers  or  the  substrate  itself. 


Figure  20.  CL  Spectra  of  1  pm  Thick  Samples  of  Each  Set  at  300K 

C.  EXPERIMENTAL  RESULTS 

The  main  objective  of  the  experiment  was  to  apply  the  transport  imaging 
technique  on  each  of  the  twelve  samples  in  order  to  determine  their  minority 
carrier  diffusion  lengths.  The  mobility  values  could  therefore  be  subsequently 
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determined  since  independent  TRPL  data  were  available.  For  all  the 
experiments  on  the  InGaP  samples,  the  SEM  was  set  with  the  experimental 
configuration  as  given  in  Table  2. 


SETTING 

VALUE 

Electron  Beam  Energy 

10  keV 

Probe  Current 

3  x  10-11  A 

Magnification 

x  500 

Table  2.  SEM  Operating  Parameters  for  InGaP  Measurements 

For  each  sample,  two  sets  of  measurements  were  taken — one  at  room 
temperature  and  the  other  at  low  temperature  (5K). 

1.  Room  Temperature  Measurements 

Figure  21  shows  the  luminescence  distributions  of  all  the  samples  from 
Set  01  to  Set  04  taken  at  room  temperature  (i.e.,  300k).  These  luminescence 
distributions  were  obtained  by  subtracting  the  ambient  light  and  detector  thermal 
noise  (captured  in  the  dark  images)  from  the  recombination  luminescence  of  the 
samples  (captured  in  the  luminescence  images).  The  distributions  were  then 
plotted  as  semi-logarithmic  graphs  of  the  normalized  luminescence  intensity,  I N , 

plotted  against  the  position,  x . 


31 


Position  [microns] 


(a)  Set  01 


(b)  Set  02 
32 


(c)  Set  03 


(d)  Set  04 


Figure  21 .  Normalized  Luminescence  Intensity  as  a  Function  of  Position  at  300K 

for  Samples  in  Sets  01  to  04 
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2.  Low  Temperature  Measurements 

Figure  22  shows  the  luminescence  distributions  of  all  the  samples  from 
Set  01  to  Set  04  taken  at  5K. 


Position  [microns] 


(a)  Set  01 


(b)  Set  02 
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(c)  Set  03 


(d)  Set  04 


Figure  22.  Normalized  Luminescence  Intensity  as  a  Function  of  Position  at  5K 

for  Samples  in  Sets  01  to  04 


35 


3.  Results  of  Diffusion  Length  Extraction 

From  the  luminescence  distribution  of  the  samples,  the  gradient  of  the 
distribution  slopes  could  be  obtained  and  the  diffusion  length  was  calculated  as 
the  inverse  of  the  slope  gradient,  using  the  1/slope  technique. 

Figure  23  illustrates  how  line-fitting  was  performed  on  the  diffusion 
distributions  for  four  of  the  samples  (two  at  300K  and  two  at  5K)  in  order  to 
determine  the  gradient  of  the  slopes.  In  general,  the  position  data  (x-axis)  used 
for  the  line-fitting  was  approximately  in  the  range  of  1 .0  pm  to  3.6  pm,  though  the 
actual  data  used  varied  for  each  set  depending  on  the  distribution.  Interestingly, 
the  distribution  of  Set  3  samples  spread  out  significantly  at  low  temperature, 
allowing  a  wider  range  of  data  to  be  used  for  the  line-fitting  (Figure  23d). 
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Figure  23.  Plots  of  Ln  (Normalized  Intensity)  as  a  Function  of  Position  for 
Gradient  Extraction  for  InGaP  Samples 
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Table  3.  shows  the  calculation  of  the  diffusion  lengths  using  the  gradient 
of  the  slopes  obtained  via  the  earlier  line-fitting  process. 


Sample 

Temperature 

(K) 

Gradient 

(Mm'1) 

Diffusion 

Length 

(pm) 

IB 

300 

-0.997 

1.00 

4A 

-1.023 

0.98 

2C 

5 

-0.892 

1.12 

3B 

-0.448 

2.23 

Table  3.  Calculation  of  Diffusion  Lengths 


The  summary  of  the  diffusion  length  values  for  all  the  twelve  samples  is 
tabulated  in  Table  4. 


Set# 

Type 

Sample 

Thickness 

(pm) 

Diffusion 
Length  @ 
300K 
(pm) 

Diffusion 
Length  @  5K 
(pm) 

1  1 

1 

p-InGaP 

(Higher  doping  -  1  x  1 017  cm"3) 

1A 

0.3 

1.11 

1.24 

IB 

0.5 

1.00 

1.18 

1C 

1.0 

0.99 

1.13 

1  1 

2 

p-InGaP 

(Lower  doping  -  5  x  1 016  cm"3) 

2A 

0.3 

1.17 

1.28 

2B 

0.5 

1.12 

1.15 

2C 

1.0 

1.06 

1.12 

1  1 

3 

p-InGaP 

(Lower  doping  -  5  x  1 016  cm"3) 

3A 

0.3 

1.60 

2.15 

3B 

0.5 

1.30 

2.23 

3C 

1.0 

1.24 

2.18 

1  1 

4 

n-InGaP 

(Doping  -  8  x  1017  cm'3) 

4A 

0.3 

0.98 

0.98 

4B 

0.5 

0.95 

1.00 

4C 

1.0 

0.94 

0.95 

Table  4.  InGaP  Sample  Diffusion  Length  Results 
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D.  ANALYSIS  OF  RESULTS 

1.  Material  Type 

The  measurements  clearly  showed  that  Set  01  to  Set  03  samples,  which 
were  p-type  materials,  had  broader  luminescence  distributions  than  Set  04 
samples,  which  were  n-type  materials.  Figure  24  shows  one  such  comparison  of 
the  distributions  for  the  set  of  samples  with  0.3  pm  InGaP  layers  at  room 
temperature.  This  indicated  that  the  p-type  materials  had  longer  diffusion  lengths 
than  the  n-type  materials.  In  p-type  materials,  the  minority  charge  carriers  were 
electrons  while  the  minority  charge  carriers  for  n-type  materials  were  holes. 
Literature  values  for  majority  carrier  mobility  have  indicated  that  the  mobility  for 
electrons  is  much  higher  than  that  of  holes  [17],  which  is  likely  to  result  in  longer 
diffusion  lengths  when  the  lifetimes  are  comparable.  Therefore,  it  is  within 
expectation  for  the  p-type  materials  to  have  longer  diffusion  lengths  than  the  n- 
type  materials. 


Figure  24.  Comparison  of  Luminescence  Distributions  of  Samples  with  0.3  pm 

InGaP  Layer  at  300K 
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2.  Doping  Level 

By  comparing  the  p-type  materials  (i.e. ,  Set  01  to  Set  03),  the 
experimental  results  showed  that  samples  in  Set  02  and  Set  03,  with  a  lower 
doping  level,  had  longer  diffusion  lengths  than  those  in  Set  01,  which  were 
higher-doped.  This  was  also  shown  by  the  broader  luminescence  distributions  of 
the  higher-doped  samples  in  Figure  24. 

In  semiconductors,  the  doping  level  has  an  influence  on  the  diffusivity  of 
the  minority  charge  carriers.  When  dopants  are  introduced  into  a  material,  they 
are  generally  ionized  at  room  temperature,  by  either  accepting  or  donating 
electrons  to  the  material.  In  either  case,  the  ions,  being  charge  centers,  act  as 
scattering  sites.  When  the  minority  charge  carriers  diffuse  through  the  material, 
they  can  be  scattered  or  deflected  from  their  original  path  by  the  Coulomb  force 
exerted  by  these  dopant  ions.  Therefore,  as  the  doping  level  increases,  which 
implies  a  higher  dopant  concentration,  the  probability  of  the  minority  charge 
carriers  encountering  scattering  or  deflection  is  now  greater.  This,  in  turn,  lowers 
the  diffusivity  of  the  charge  carriers,  resulting  in  shorter  diffusion  lengths. 

3.  Thickness 

From  the  graphs  (Figure  21  for  room  temperature  measurements  and 
Figure  22  for  low  temperature  measurements),  it  could  be  seen  that  in  general, 
the  distribution  narrowed  with  increasing  thickness  for  each  set.  The  difference 
in  the  distribution  width  was  very  clearly  demonstrated  in  the  case  of  Set  02  and 
Set  03  samples,  though  it  was  not  so  obvious  in  the  case  of  Set  01  and  Set  04 
samples.  Since  a  narrower  distribution  has  steeper  slopes,  the  slope  gradient 
will  have  a  higher  value  when  compared  with  that  of  broader  distributions.  This 
implies  that  the  diffusion  length,  which  was  calculated  as  the  inverse  of  the  slope 
gradient,  would  be  shortest  for  the  sample  with  the  highest  thickness  within  each 
set. 

A  possible  explanation  for  this  is  that  as  the  thickness  increases,  the 
minority  charge  carriers  are  less  confined  to  ID  or  2D  diffusion  and  start  to 
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experience  diffusion  in  the  depth  too,  resulting  in  3D  diffusion.  As  such,  the 
charge  carriers  diffuse  in  a  volumetric  manner  and  the  apparent  diffusion  length 
decreases,  since  the  2D  model  was  applied  in  all  cases. 

4.  Temperature 

A  comparison  of  the  luminescence  distributions  of  sample  1C  at  various 
temperatures  (i.e.,  from  300K  to  5K)  is  given  in  Figure  25.  From  the  graph,  it 
could  be  observed  that  as  temperature  dropped,  the  distribution  started  to  widen, 
thereby  giving  longer  diffusion  lengths  for  the  same  sample.  This  observation 
was  the  same  for  the  other  sample  sets. 

The  diffusivity  of  minority  charge  carriers  is  also  influenced  by 
temperature.  Due  to  thermal  energy,  the  atoms  within  a  lattice  are  not  at  rest, 
but  constantly  vibrating  about  their  fixed  positions.  These  quantized  vibrational 
states  are  known  as  phonons.  As  the  charge  carriers  travel  through  the  material, 
they  will  eventually  interact  with  such  a  phonon  and  be  deflected  from  their 
original  paths.  Thus,  the  vibration  causes  scattering  of  charge  carriers.  As 
temperature  increases,  the  atoms  gain  more  thermal  energy  and  the  population 
of  phonons  increases  as  well.  This  causes  more  scattering  of  the  charge 
carriers,  leading  to  lower  carrier  mobility  at  higher  temperature. 

The  actual  change  in  mobility  is  also  affected  by  the  fact  that  these 
materials  were  ternary  alloys.  In  this  case,  alloy  disorder  scattering  plays  a 
significant  role,  weakening  the  effect  of  temperature,  compared  to  a  binary 
material,  such  as  GaAs  or  Indium  Phosphide  (InP). 
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Figure  25.  Comparison  of  Sample  1C  Luminescence  Distributions  at  Various 

Temperatures 


E.  RESULTS  OF  MOBILITY  VALUE  CALCULATION 


The  next  objective,  having  obtained  the  diffusion  lengths,  was  to 
determine  the  mobility  values  of  the  minority  charge  carriers  in  the  samples.  The 
sample  providers  had  earlier  acquired  the  minority  charge  carrier  lifetime  data  of 
their  samples  by  means  of  TRPL.  Using  the  diffusion  length  data,  together  with 
the  lifetime  data  set,  the  minority  charge  carrier  mobility,  // ,  could  be  easily 
calculated  using  the  following  equation: 


M  = 


ikT 


Eq.  (13) 


where  z  is  the  diffusion  length,  e  is  the  electron  charge,  r  is  the  minority 


charge  carrier  lifetime,  k  is  the  Boltzmann’s  constant  and  T  is  the  temperature. 
Note  that  this  is  just  a  rearrangement  of  Equation  3.  It  should  be  highlighted  that 
the  lifetime  data  were  obtained  at  room  temperature  settings  and,  therefore, 
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could  not  be  used  to  evaluate  the  mobility  values  of  the  samples  at  low 
temperature.  The  results — for  the  room  temperature  measurements  only — are 
summarized  in  Table  5. 


Set# 

Type 

Sample 

Thickness  (pm) 

Diffusion 
Length  @  300K 
(pm) 

Lifetime  @ 
300K 
(ns) 

Mobility 

@ 

300K 

(crrPV's 

f) 

1 

p-InGaP 
(higher  doping  - 
1  x  1 017  crrr3) 

1A 

0.3 

1.11 

1.13 

421 

IB 

0.5 

1.00 

1.03 

375 

1C 

1.0 

0.99 

1.32 

287 

2 

p-InGaP 

(lower  doping  -  5 
x  1016  cm-3) 

2A 

0.3 

1.17 

1.07 

494 

2B 

0.5 

1.12 

0.83 

584 

2C 

1.0 

1.06 

0.95 

457 

3 

p-InGaP 

(lower  doping  -  5 
x  1016  cm-3) 

3A 

0.3 

1.60 

1.36 

727 

3B 

0.5 

1.30 

1.37 

477 

3C 

1.0 

1.24 

1.55 

383 

4 

n-InGaP 
(higher  doping  - 
1  xIO17  cm-3) 

4A 

0.3 

0.98 

228 

163 

4B 

0.5 

0.95 

2.07 

168 

4C 

1.0 

0.94 

2.11 

162 

Table  5.  InGaP  Minority  Charge  Carrier  Mobility  Results 

F  OBSERVATIONS 

One  observation  of  the  TRPL  data  set  was  that  the  minority  carrier 
lifetimes  did  not  follow  a  specific  trend.  Normally,  one  would  expect  the  lifetime 
values  to  increase  with  increasing  thickness  for  samples  within  each  set  of  the 
same  material  type  and  doping  level  due  to  improvement  in  material  quality  and  a 
smaller  role  in  surface  recombination.  However,  this  was  not  really  observed 
and,  instead,  there  were  many  cases  where  the  lifetimes  were  shorter  for  the 
thicker  samples.  This  showed  that  the  fabrication  process  employed  to 
manufacture  these  samples  might  not  have  been  mature  yet,  resulting  in 
inconsistencies  in  the  material  uniformity  or  defects,  which  could  affect  the 
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transport  behavior  of  the  charge  carriers.  Another  possibility  was  that  surface 
recombination  in  these  samples  was  very  low  and  the  variations  in  lifetime  data 
represented  the  experimental  error. 

Another  critical  observation  from  the  results  was  that  the  hole  mobility 
values  extracted  from  the  simple  slope  analysis  appeared  to  be  too  high. 
Though  literature  data  on  minority  charge  carrier  mobility  is  scarce,  there  is  good 
indication  that  the  norm  of  such  mobility  is  in  the  range  of  100  to  1000  cm2V"1s"1 
for  electrons  and  10  to  100  cmVs'1  for  holes  [17].  While  the  electron  carriers 
presented  mobility  values  of  300  to  700  cm2V"1s"1,  which  were  within  the 
expected  range,  the  experimental  results  of  about  160  cmVs'1  for  holes  were 
higher  than  expected. 

One  possibility  for  this  phenomenon  was  that  the  actual  diffusion  lengths 
of  the  samples  were  much  shorter  than  those  that  were  derived  from  the  1  /slope 
technique  of  the  luminescence  distribution.  Due  to  the  short  diffusion  lengths, 
the  luminescence  distributions  seen  in  the  graphs  were  most  likely  a  convolution 
of  the  actual  diffusion  behavior  with  the  finite  size  of  the  SEM’s  electron  beam. 
Therefore,  the  diffusion  lengths  of  the  samples  could  not  be  effectively 
determined  by  using  the  1/slope  technique  of  this  convolved  luminescence,  since 
it  was  also  affected  by  the  size  of  the  generation  region  for  the  charge  carriers. 

In  earlier  work  on  transport  imaging,  the  samples  all  had  considerably  long 
diffusion  lengths,  such  that  the  recombination  luminescence  occurred  even  when 
“far  away”  from  the  electron  beam.  This  allowed  the  1/slope  technique  to  remain 
effective  as  the  luminescence  distribution  was  dominated  by  diffusion  without 
major  interference  from  the  finite  generation  region.  For  these  samples, 
however,  the  1/slope  technique  failed.  To  determine  the  real  diffusion  lengths, 
there  is  now  a  need  to  separate  the  resulting  luminescence  into  its  individual 
components — the  diffusion  behavior  in  the  samples  and  the  incident  electron 
beam  interaction.  This  will  be  discussed  in  the  next  chapter. 
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V.  DECONVOLUTION  MODEL  FOR  OVERCOMING  THE  FINITE 
VOLUME  GENERATION  LIMITATION 

A.  EXPERIMENTAL  LIMITATION 

The  results  of  the  charge  carrier  mobility  values  presented  in  the  previous 
chapter  were  found  to  be  significantly  higher  than  would  be  expected  based  on 
literature  values  for  majority  carrier  mobility  in  InGaP,  particularly  for  the  minority 
holes  in  the  n-type  samples.  One  hypothesis  was  that  the  minority  diffusion 
lengths  of  the  samples  were  sufficiently  short  (0.5  pm  to  2.0  pm)  to  be  of  the 
same  order  of  magnitude  as  the  finite  charge  carrier  volume  generation  region 
due  to  the  electron  beam.  As  such,  when  the  transport  imaging  technique  was 
applied,  what  was  captured  by  the  CCD  camera  was  the  convolved  product  of 
the  actual  recombination  luminescence  from  the  samples  with  the  finite  charge 
carrier  volume  generation.  Therefore,  when  the  charge  carrier  mobility  was 
calculated  without  including  the  effect  of  the  finite  size  of  the  source  beam,  the 
results  deviated  significantly. 

The  volume  generation  region  can  be  illustrated  using  a  Monte  Carlo 
simulation.  CASINO,  a  simulation  program  based  on  the  Monte  Carlo  algorithm, 
was  used  to  simulate  the  effect  of  transport  imaging  using  a  10  keV  beam  source 
on  a  0.3  pm  thick  InGaP  layer  on  GaAs  (Figure  26).  From  the  graph,  the 
existence  of  the  generation  region  was  clearly  depicted  and  the  region  even 
extended  to  a  depth  beyond  the  InGaP  material.  Laterally,  the  generation  region 
extended  to  a  width  of  approximately  0.5  pm.  The  consequence  is  that  if  the 
diffusion  length  of  the  samples  were  short,  the  recombination  luminescence 
would  become  coupled  with  the  effect  of  the  generation  volume.  This  is  a 
limitation  that  must  be  addressed  in  the  model  for  short  diffusion  length  materials. 
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Figure  26.  Monte  Carlo  Simulation  Showing  Finite  Generation  Volume.  The 
Image  Shows  the  Energy  Distribution  in  Both  Lateral  and  Depth  Dimensions. 


Figure  27  shows  the  effect  of  convolution  between  the  generation  volume 
region  and  the  recombination  luminescence.  Without  taking  this  into  account,  the 
diffusion  length  measured  by  the  transport  imaging  will  not  be  a  correct 
interpretation  of  the  actual  minority  charge  carrier  diffusion  behavior. 


Figure  27.  Convolution  of  the  Generation  Volume  with  the  Recombination 

Luminescence 
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To  overcome  this  limitation  and  obtain  the  true  mobility  values,  there  was 
a  need  to  deconvolve  the  resultant  luminescence  back  into  its  original 
components  so  that  they  could  be  evaluated  independently. 


B.  DECONVOLUTION  MODEL 


The  ideal  method  was  to  perform  a  direct  deconvolution  process  on  the 
results  to  retrieve  the  actual  recombination  luminescence.  To  do  this,  there  was 
a  need  to  model  the  two  components — both  the  generation  volume  and  the 
diffusion  behavior — in  mathematical  form. 

In  this  first  attempt,  the  finite  charge  carrier  generation  volume  could  be 
approximated  as  a  region  that  followed  a  Lorentz  distribution  perpendicular  to  the 
line  source  and  could  be  represented  as: 


L{x)  = 


Y 

(x-x0)2-y2 


Eq.  (14) 


where  x  is  the  horizontal  position,  x0  is  the  center  of  the  distribution,  and  y  is 

the  Lorentz  parameter  characterizing  the  width  of  the  generation  region.  Since 
the  measurements  were  taken  from  the  center,  x0  was  set  to  0  and  Equation  14 


reduced  to: 


L{x) 


Y 


2  2 

x  -y 


Eq.  (15) 


From  the  information  presented  in  Chapter  III,  the  diffusion  distribution 
from  a  point  source  could  be  described  by: 


IJx)=e 


w 


Eq.  (16) 


where  x  is  the  horizontal  position  and  Ldiff  is  the  diffusion  length. 

The  convolved  result,  R ,  of  the  carrier  diffusion  from  the  charge  carrier 
generation  volume  region  could  now  be  represented  as: 


Y 


e+Y2 


-• e 


x-z\ 

Ldiff 


Eq.  (17) 


where  C,  is  the  integration  variable  representing  position  in  the  diffusion  axis. 
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It  should  be  noted  that  while  this  would  be  the  most  direct  method  to 
independently  determine  Ldifr  and  y ,  it  would  require  a  two-parameter  fit  of  the 

integrated  term.  Due  to  the  difficulty  of  performing  direct  deconvolution,  a 
proposed  two-stage  deconvolution  process  would  be  employed  instead.  This 
two-stage  process  involved  the  using  of  known  data  to  first  approximate  the 
value  of  the  y  parameter,  and  thereafter  to  calculate  the  diffusion  lengths  based 
on  this  approximated  y  parameter.  This  process  was  proposed  as  a  means  of 
testing  the  viability  of  the  deconvolution  model  while  the  actual  multi-parameter 
fitting  could  be  carried  out  in  the  future  if  the  approach  proved  to  be  feasible. 

C.  DECONVOLUTION  PROCESS 

1 .  Stage  I 

The  first  stage  sought  to  approximate  an  empirical  description  of  the 
generation  region  due  to  the  electron  beam  (i.e.,  the  y  parameter  of  the  Lorentz 
distribution),  and  then  determine  the  diffusion  length  based  on  this  starting  point. 
To  achieve  this,  the  mobility  value  for  one  sample  set  had  to  be  assumed  in  order 
to  estimate  the  value  of  y  parameter.  To  reduce  the  uncertainty  of  the  assumed 
mobility  value,  the  sample  with  the  supposed  shortest  diffusion  length  would  be 
chosen.  This  would  be  the  Set  04  samples. 

In  terms  of  minority  charge  carrier  mobility  values,  Set  04  samples  should 
theoretically  be  the  lowest  among  all  four  sample  sets.  The  first  reason  was  that 
Set  04  samples  were  n-type  materials,  thus  having  holes  as  their  minority  charge 
carriers  as  opposed  to  electrons  for  the  others.  From  Hall  effect  measurements 
[17],  it  has  been  shown  that  holes  in  InGaP  have  lower  mobility  than  electrons, 
thus  contributing  to  shorter  diffusion  lengths.  The  second  reason  was  that  it  was 
also  the  sample  set  with  the  highest  doping  level.  As  discussed  earlier  in 
Chapter  IV,  a  higher  doping  level  generally  reduces  the  diffusion  length  due  to 
lower  mobility  and,  in  many  cases,  shorter  lifetimes.  Therefore,  being  an  n-type 
material  and  being  the  heaviest  doped  sample  set,  Set  04  samples  were 
expected  to  have  the  shortest  minority  carrier  diffusion  lengths. 
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Based  on  the  Hall  effect  on  the  mobility  of  holes  as  majority  charge 
carriers  [18],  the  mobility  value  of  holes  as  minority  charge  carrier  was  assumed 
to  be  45  cm2V1s"1.  With  the  measured  lifetime  data  (r  =  2.28  ns),  the  diffusion 
length  was  calculated  to  be  0.52  pm. 

Using  the  approximated  diffusion  length  value  for  z,  in  Equation  17, 
different  values  of  y  were  substituted  and  the  predicted  intensity  distributions  of 
R  were  calculated  (Figure  28).  From  the  R  distributions,  linear  regression  was 
performed  to  obtain  the  gradient  of  the  slope  (Figure  29).  The  range  of  data 
used  for  the  line-fitting  was  in  the  range  of  0.4  to  0.8  for  normalized  intensity  as 
this  was  the  part  where  the  curve  stabilized  and  presented  a  relatively  straight 
profile.  The  convolved  diffusion  lengths  of  these  R  results  were  then  compared 
against  the  actual  experimental  result  for  Sample  4A  to  identify  the  case  that 
gave  the  closest  fit.  The  y  parameter  value  was  then  selected  to  be  the  y  value 
that  gave  this  closest  fitting  result.  After  iterations  of  y  value  substitution  and 
comparison,  the  y  parameter  value  was  determined  to  be  0.48  pm  (Table  6). 
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Figure  28.  Results  for  Different  Values  of  the  Lorentz  Parameter 
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Figure  29.  Line-fitting  for  Extraction  of  Slope  Gradient  for  Stage  I 


Lorentz 

Parameter 

(pm) 

Slope  Gradient 
(Mm1) 

1/Slope  Gradient 
(pm) 

Experimental 
Diffusion  Length 
(pm) 

0.08 

-1.725 

0.58 

0.28 

-1.318 

0.76 

0.48 

-1.021 

0.98 

0.98 

0.68 

-0.814 

1.23 

0.88 

-0.664 

1.51 

Table  6.  1/Slope  Gradient  Results  for  Different  Values  of  the  Lorentz 

Parameter 


2.  Stage  II 

Having  determined  the  best  estimated  value  of  the  y  parameter,  the 
deconvolution  model  now  became  a  one-parameter  fit  model — that  of  z  . 
Using  the  value  of  0.48  pm  for  the  y  parameter,  the  Lorentz  distribution 
representing  the  generation  volume  was  convolved  with  the  exponential 
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distribution  of  the  recombination  luminescence  with  different  values  of  Ldiff  for 
each  sample.  Similarly,  the  z  value  for  each  sample  was  determined  based 

on  the  value  that  gave  the  closest  match  with  the  actual  experimental  results. 
Figure  30  illustrates  the  different  resulting  distributions  based  on  substituting 
different  values  of  z  for  Sample  3A  and  Figure  31  shows  the  line-fitting 

process  for  the  extraction  of  the  slope  gradient  of  the  distributions.  The  diffusion 
lengths  of  the  resulting  distributions  for  Sample  3A  were  then  calculated  and 
tabulated  in  Table  7.  The  case  of  Ldiff ,  with  a  value  of  1 .22  pm,  gave  the  closest 

match  of  the  convolved  diffusion  length  with  that  of  experimental  results. 
Therefore,  the  Ldiff  value  for  Sample  3A  is  determined  to  be  1 .22  pm. 
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Figure  30.  Results  for  Different  Values  of  the  Diffusion  Length  for  Sample  3A 
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Figure  31 .  Line-fitting  for  Extraction  of  Slope  Gradient  for  Stage  II 


Actual 

Diffusion 

Length, 

(Mm) 

Slope  Gradient 
(Mm'1) 

1/Slope  Gradient 
(pm) 

Experimental 
Diffusion  Length 

(pm) 

1.02 

-0.706 

1.42 

1.22 

-0.625 

1.60 

1.60 

1.44 

-0.558 

1.79 

Table  7.  1 /Slope  Gradient  Results  for  Different  Values  of  the  Diffusion  Length 


D.  RESULTS  AFTER  DECONVOLUTION 

The  application  of  the  deconvolution  model  was  able  to  separate  the 
effects  of  the  finite  charge  carrier  generation  volume  from  that  of  the  diffusion 
behavior,  thereby  obtaining  more  accurate  diffusion  lengths  for  each  of  the 
samples.  The  recalculated  mobility  of  the  minority  charge  carriers  based  on  the 
deconvolved  diffusion  lengths  were  reduced  in  value  and  were  now  order-of- 
magnitude  consistent  with  the  limited  expected  values  based  on  other  techniques 
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in  the  literature.  The  deconvolved  diffusion  lengths  and  the  recalculated  minority 
charge  carrier  mobility  values  are  summarized  in  Table  8. 


Set# 

Type 

Thickness 

(pm) 

Lifetime 
@  300K 
(ns) 

Uncorrected 
Diffusion 
Length  @  300K 
(pm) 

Deconvolved 
Diffusion 
Length  @  300K 
(pm) 

Uncorrected 
Mobility  @ 
300K 

(cmW) 

Deconvolved 
Mobility  @ 
300K 

(cm2V*1s*1) 

1 

p-InGaP 

(higher  doping  - 
1  x  1017  cm*3) 

0  3 

1.13 

1.11 

0.67 

421 

154 

0.5 

1.03 

1.00 

0.56 

375 

118 

1.0 

1.32 

0  99 

0.54 

287 

85 

2 

p-InGaP 

(lower  doping-  5 
x  10"  cm-3) 

0  3 

1.07 

1.17 

0  75 

494 

203 

0.5 

0.83 

1.12 

0  68 

584 

215 

1.0 

0.95 

1.06 

0.62 

457 

156 

3 

p-InGaP 
(lower  doping  - 
5  x  10"  cnr3) 

0.3 

1.36 

1.60 

1.22 

727 

423 

0.5 

1.37 

1  30 

0  90 

477 

228 

1.0 

1.55 

1.24 

0.83 

383 

172 

Table  8.  InGaP  Minority  Charge  Carrier  Mobility  Results  After  Deconvolution 

E.  SENSITIVITY  TEST 

A  sensitivity  analysis  was  also  performed  to  determine  the  degree  of 
variation  in  the  results  due  to  different  assumptions  for  the  starting  value  of  the 
hole  mobility  for  the  deconvolution  process.  Earlier,  the  assumed  value  for  the 
hole  mobility  was  taken  to  be  45  cm2V-1s-1.  For  the  sensitivity  analysis,  two 
new  hole  mobility  values  were  used  (one  was  double  of  the  original  value  while 
the  other  was  half  the  value)  and  the  deconvolution  process  re-run  to  determine 
the  new  z  for  the  samples.  The  results  based  on  the  diffusion  of  Sample  3A 

are  given  in  Table  9. 
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Assumed 

Hole  Mobility  @ 
300K 

(cmW) 

Resultant 

Lorentz 

Parameter 

(Mm1) 

Diffusion  Length 
(pm) 

Electron 
Mobility  @ 
300K 

(cmW) 

22.5 

0.56 

1.14 

369 

45.0 

0.48 

1.22 

423 

90.0 

0.31 

1.37 

533 

Table  9.  Sensitivity  Testing  of  Sample  3A 

From  the  sensitivity  analysis  results,  it  was  observed  that  when  the 
assumed  hole  mobility  value  was  halved  (i.e.,  22.5  crr^vV1),  the  electron 
mobility  value  was  decreased  by  13%.  Similarly,  when  the  hole  mobility  was 
assumed  to  be  double  the  base  value  (i.e.,  90  cmV's'1),  there  was  an  increase 
of  26%  to  the  mobility  of  the  electrons.  From  a  statistical  point  of  view,  the 
electron  mobility  could  be  considered  rather  sensitive  to  the  starting  value  of  the 
hole  mobility  used  for  the  deconvolution  process  and  could  potentially  vary  by 
13%  to  26%.  However,  from  an  absolute  value  point  of  view,  the  results  were 
deemed  to  be  stable  as  the  numerical  change  is  small  when  compared  to  the 
mobility  values  of  majority  carriers,  which  are  in  the  range  of  1000  cm2V"1s"1  or 
higher,  and  were  order-of-magnitude  consistent  with  the  expected  values  from 
literature  data. 

F.  OBSERVATIONS 

The  range  of  the  majority  charge  carrier  mobility  (for  both  electrons  and 
holes)  based  on  the  Hall  effect  was  also  provided  by  the  sample  providers  and  is 
given  in  Table  10.  A  comparison  between  Table  8  and  Table  10  showed  that  the 
mobility  of  the  electrons  as  minority  carriers  was  considerably  lower  than  that  of 
the  electrons  as  majority  carriers.  This  suggests  that  minority  charge  carrier 
mobility  is  lower  than  majority  charge  carrier  mobility  for  a  given  carrier  type  in  a 
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single  material.  The  finding  is  consistent  with  theoretical  work  in  GaAs,  which 
predicts  that  minority  carrier  mobilities  are  lower  by  a  factor  of  0.6  and  0.9, 
depending  on  compensation  [19]. 


Material 

Majority  Charge 
Carrier  Type 

Charge  Carrier 
Concentration 

(cm'3) 

Hall  Mobility 
(cmW) 

1  x  10'IB 

800 

Electron 

2  x  101B 

650 

5  x  101B 

500 

InGaP 

7  x  10  ia 

450 

2x10™ 

45 

Hole 

3x10™ 

40 

8  x  101B 

27 

Table  10.  InGaP  Majority  Charge  Carrier  Mobility  Based  on  Hall  Effect  (From 

[18]) 


Based  on  the  results  of  the  sensitivity  analysis,  it  was  also  shown  that  the 
deconvolution  model  was  fairly  successful  in  eliminating  the  effect  of  the 
generation  volume  from  the  recombination  luminescence  to  allow  the  extraction 
of  the  actual  diffusion  length  and,  subsequently,  the  minority  charge  carrier 
mobility.  Although  the  method  presented  in  this  thesis  was  a  first  attempt  at  the 
deconvolution  process,  it  has  demonstrated  that  the  model  is  a  feasible  one  and 
could  prove  more  effective  with  further  research. 
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VI.  TRANSPORT  IMAGING  OF  CIGS 


A.  CIGS  TECHNOLOGY 

Copper  Indium  Gallium  (di)Selenide,  also  known  as  CulnGaSe2  or  CIGS, 
represents  an  alternate  semiconductor  material  choice  that  can  be  used  for  solar 
cell  manufacturing.  CIGS  is  a  l-lll-VI2  semiconductor  that  is  made  of  the 
elements  copper,  indium,  gallium  and  selenium,  as  the  name  will  have  implied. 
The  elements  are  tetrahedrally  bonded  in  a  chalcopyrite  crystal  structure  to  form 
the  CIGS  compound.  It  is  a  direct  band-gap  material  with  a  band-gap  range  of 
1.04  eV  to  1.67  eV  depending  on  the  composition  of  indium  and  gallium  mixture 
in  the  compound  [20], 

CIGS  suffers  from  low  energy-conversion  efficiency  when  compared  to 
other  solar  cell  materials.  Currently,  the  highest  reported  conversion  efficiency 
stands  at  20.3%,  while  the  average  efficiency  attained  for  normal  line  production 
standard  is  at  10  to  12%.  Alongside  single-junction  silicon-based  solar  cells  with 
peak  efficiency  of  27.6%  and  an  average  of  16  to  17%  efficiency,  or  more 
advanced  multi-junction  solar  cells  of  30  to  40%  efficiency  [21],  CIGS  pales  in 
comparison. 

CIGS  does,  however,  have  a  few  unique  characteristics  that  are 
advantageous  over  other  material  types,  and  these  are  the  reasons  why  solar 
cell  producers  are  interested  in  research  and  commercial  development  of  CIGS 
material. 

One  such  advantage  of  CIGS  over  silicon-based  or  other  material  type 
semiconductors  is  that  it  has  a  very  high  absorption  coefficient  (105  cm'1  for  1.5 
eV  or  higher  energies)  as  compared  to  the  rest  [22],  Up  to  99%  of  light  incident 
of  it  is  absorbed  in  the  first  1  pm  or  so  [23],  The  implication  is  that  only  a  thin  film 
of  CIGS  is  needed  to  absorb  the  same  amount  of  sunlight  compared  to  other 
semi-conductor  materials  with  lower  absorption  coefficients. 

Another  advantage  of  CIGS  is  that  it  can  be  deposited  on  flexible 
substrates,  such  as  polymide  or  metal  foils  [24],  Coupled  with  its  high  absorbing 
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characteristic,  CIGS  is  an  ideal  choice  for  the  manufacturing  of  thin-film  flexible 
solar  cells,  which  are  relatively  lightweight  and  can  be  easily  portable.  Unlike 
crystalline  silicon,  which  produces  rigid  and  delicate  solar  panels  that  require 
careful  handling,  these  flexible  and  portable  CIGS  solar  cells  offer  a  myriad  of 
product  opportunities  in  both  the  commercial  and  military  sectors.  One  possible 
application  in  the  military  domain  is  a  CIGS  solar  cell  mat  that  can  be  rolled  up 
and  carried  by  the  average  soldier  during  movement  and  can  be  easily  set  up  to 
power  devices,  such  as  small-scale  unmanned  aerial  vehicles  (UAVs)  or 
communication  sets,  during  operations. 

One  deciding  factor  of  choosing  CIGS  though,  will  be  in  terms  of  cost. 
CIGS  is  a  relatively  cheaper  material  to  manufacture  and  install  compared  to 
mainstream  solar  cell  materials,  such  as  silicon  or  InGaP.  This  is  due  to  its 
capacity  of  operating  with  polycrystalline,  as  opposed  to  higher-cost  single¬ 
crystal  materials,  thus  making  it  a  lower-cost  material.  The  CIGS  technology  of 
today  is  far  from  mature  but  with  further  research,  the  production  cost  of  CIGS- 
based  solar  cells  can  be  driven  down  even  lower.  Together  with  its  other 
interesting  characteristics  of  high  absorbency  and  substrate  flexibility,  CIGS  solar 
cells  have  the  potential  to  grow  and  become  a  huge  competitor  in  the  solar  cell 
industry. 

B.  SAMPLES 

Through  collaboration  with  researchers  in  the  solar  cell  field,  three 
samples  of  CIGS  solar  cells  were  obtained  for  studying  in  this  thesis  work. 
These  samples  will  be  labeled  as  Sample  Cl,  C2,  and  C3,  respectively. 

Transport  imaging  has  been  performed  on  these  three  samples  with  the 
objective  of  measuring  their  diffusion  lengths  and  finding  their  minority  charge 
carrier  mobility  value  in  order  to  understand  their  transport  behavior  better.  As 
CIGS  are  p-type  materials,  the  focus  is  the  electron  carrier  mobility. 
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C.  SEM  SETTINGS 

For  all  the  experiments  on  the  CIGS  samples,  the  SEM  was  set  with  the 
experimental  configuration  as  given  in  Table  1 1 : 


SETTING 

VALUE 

Electron  Beam  Energy 

20  keV 

Probe  Current 

3  x  10'10  A 

Magnification 

x  1500 

Table  1 1 .  SEM  Operating  Parameters  for  CIGS  Measurements 

D.  INITIAL  TRANSPORT  IMAGING  OF  CIGS 

The  luminescence  images  of  the  CIGS  samples  using  the  transport 
imaging  technique  are  given  in  Figure  32.  In  contrast  to  the  luminescence  seen 
in  the  earlier  InGaP  samples,  the  luminescence  observed  of  the  CIGS  had  less 
uniform  distributions.  The  distorted  luminescence  pattern  was  either  a  reflection 
of  material  non-uniformity  or  a  result  of  light-scattering  caused  by  the  rough  and 
uneven  surface  of  the  material. 


Figure  32.  Live-scan  Luminescence  Images  of  CIGS  Samples  [256  pixels  (102 

pm)  x  172  pixels  (68  pm)  each] 
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In  order  to  obtain  a  more  accurate  interpretation  of  the  luminescence 
distribution,  each  of  the  luminescence  images  had  to  be  separated  into  smaller 
parts  (Figure  33).  This  allowed  the  luminescence  to  be  localized  to  a  region  for 
study  and  the  various  regions  could  later  be  compared. 


Figure  33.  Localized  Luminescence  Distributions 

Figure  34  shows  a  comparison  of  the  obtained  luminescence  distributions 
of  each  of  the  samples.  From  the  graph,  it  could  be  observed  that  Sample  Cl 
had  the  widest  distribution  while  Sample  C2  had  the  narrowest. 
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Figure  34.  Comparison  of  CIGS  Luminescence  Distributions 

The  line-fitting  process  is  illustrated  in  Figure  35,  where  linear  regression 
was  performed  in  order  to  extract  the  gradient  of  the  slopes.  For  the  CIGS 
material,  a  longer  range  of  data  could  be  used  for  the  line-fitting  due  to  its  wider 
distributions  as  compared  to  those  of  the  InGaP  material.  In  general,  the  position 
data  (x-axis)  used  for  the  line-fitting  was  approximately  in  the  range  of  2  pm  to  8 
pm.  Using  the  gradient,  the  diffusion  lengths  of  the  samples  were  then 
calculated  and  tabulated  in  Table  12. 
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Figure  35.  Plot  of  Ln  (Normalized  Intensity)  as  a  Function  of  Position  for 
Gradient  Extraction  for  Sample  C3 


Samples 

Shortest 
Diffusion  Length 
(pm) 

Longest 

Diffusion  Length 
(pm) 

Average 
Diffusion  Length 
(pm) 

Cl 

3.17 

3.75 

3.49 

C2 

2.29 

2.76 

2.46 

C3 

3.38 

3.38 

3.38 

Table  12.  CIGS  Sample  Diffusion  Length  Results 

Independent  TRPL  lifetime  data  were  also  available  from  the  sample 
providers,  allowing  the  calculation  of  the  minority  charge  carrier  mobility,  which  is 
shown  in  Table  13.  Unlike  the  InGaP  samples,  there  was  no  need  for 
deconvolution  of  the  CIGS  diffusion  length.  This  was  due  to  the  fact  that  the 
diffusion  lengths  were  longer  for  the  CIGS  (i.e.,  2  pm  and  above),  and  thus  were 
out  of  the  affected  region  of  the  finite  generation  volume. 
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Samples 

Average 
Diffusion  Length 
(pm) 

Lifetime 

(ns) 

Mobility 

(cmW) 

Cl 

3.49 

4 

1180 

C2 

2.46 

6 

390 

C3 

3.38 

6 

736 

Table  13.  CIGS  Minority  Charge  Carrier  Mobility  Results 


E.  OBSERVATIONS 

From  this  experiment,  preliminary  values  for  the  minority  charge  carrier 
(electron)  mobility  of  the  CIGS  material  were  established.  While  the  mobility 
values  of  samples  C2  and  C3  were  within  expectation;  however,  the  mobility 
value  of  sample  Cl  seemed  to  be  slighter  higher  than  the  norm. 

From  the  collected  luminescence  images,  it  was  clear  that  the  surface 
roughness  played  a  part  in  scattering  the  light  source  and  distorted  the  eventual 
charge  carrier  generation  process.  This  resulted  in  the  uneven  recombination 
luminescence  that  was  observed.  It  was  believed  that  in  addition  to  the 
recombination  luminescence,  the  CCD  camera  would  have  also  received  some 
of  the  back-scattered  light.  These  coupled  signals  would  certainly  affect  the 
accuracy  of  the  transport  imaging  technique  and  hence  influence  the  subsequent 
determination  of  the  diffusion  length.  Therefore,  these  mobility  values  would  just 
represent  a  first  approximation  of  the  CIGS  transport  behavior  and  should  not  be 
taken  as  absolute  values. 
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VII.  CONCLUSION 


A.  SUMMARY 

Solar  energy  represents  a  unique  solution  to  the  world’s  increasing  energy 
demands.  To  that  end,  higher-efficiency  solar  cells  are  required  in  order  to  be 
able  to  fully  harness  the  energy  from  the  sun.  An  understanding  of  the  transport 
behavior  of  solar  cell  materials  enables  manufacturers  to  characterize  the 
materials  better  and  hence  build  more  effective  solar  cell  devices.  The  minority 
charge  carrier  diffusion  length  and  mobility  are  two  important  transport 
parameters  for  such  material  characterization. 

In  this  thesis,  work  was  done  to  determine  the  minority  charge  carrier 
diffusion  length  of  two  materials,  InGaP  and  CIGS,  which  represent  two  very 
different  solar  cell  technologies  in  the  market  today.  Transport  imaging,  a  non¬ 
destructive  and  contact-free,  optical  imaging  technique  was  applied  to  the  two 
material  samples  in  order  to  extract  the  diffusion  lengths.  Together  with 
independent  sets  of  TRPL  lifetime  data,  the  mobility  of  the  minority  charge 
carriers  were  subsequently  determined. 

For  the  InGaP  sample  set,  the  diffusion  length  was  found  to  be  in  the 
range  of  0.5  to  1.2  pm,  while  the  mobility  value  of  the  minority  charge  carriers 
was  found  to  be  in  the  range  of  80  to  400  crr^V'V1.  The  mobility  was  lower  in 
value  as  compared  to  that  of  the  majority  charge  carriers  which  was  obtained 
through  the  Hall  effect.  This  result  is  consistent  with  the  existing  literature  data  of 
GaAs  [19],  In  addition,  the  experiments  also  showed  how  the  diffusion  length 
and  mobility  values  were  affected  as  a  result  of  different  material  type,  doping 
level,  sample  thickness,  and  the  ambient  temperature. 

A  limitation  in  the  optical  technique  was  encountered  while  experimenting 
on  the  InGaP  samples.  As  the  diffusion  lengths  of  the  InGaP  samples  were 
significantly  short,  the  charge  carrier  recombination  process  was  within  the 
region  of  the  finite  generation  volume.  This  caused  the  coupling  of  the  two 
signals,  which  affected  the  initial  measurements  of  the  transport  imaging.  To 
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solve  the  problem,  a  deconvolution  process,  which  modeled  the  convolution  of 
the  recombination  luminescence  with  the  generation  volume,  was  proposed  and 
presented.  The  process  aimed  to  separate  the  two  coupled  components  into 
their  individual  parts  so  that  the  recombination  luminescence  can  be  studied 
independently.  Following  the  deconvolution  process,  both  the  diffusion  lengths 
and  minority  carrier  mobility  values  were  lowered,  allowing  a  more  accurate 
measurement  of  the  transport  parameters. 

For  the  CIGS  samples,  the  diffusion  length  was  found  to  be  in  the  range  of 
2  to  4  pm,  while  the  mobility  value  of  the  minority  charge  carriers  was  found  to  be 
in  the  range  of  300  to  1400  cm2V"1s"1.  The  challenge  of  applying  transport 
imaging  on  the  CIGS  was  that  the  rough  surface  of  the  samples  tends  to  scatter 
the  light  source  and  distorts  the  final  luminescence  image  collected  by  the  CCD 
camera.  As  a  result,  the  extracted  diffusion  lengths  based  on  the  distorted 
luminescence  represented  only  an  approximation  rather  than  the  actual  value  of 
the  transport  behavior  in  the  material. 

B  FUTURE  RESEARCH 

The  deconvolution  model  presented  in  this  thesis  has  the  objective  of 
separating  two  convolved  signals  into  individual  components  so  that  the 
luminescence  resulting  from  the  material  can  be  studied  independently  for  a 
more  accurate  interpretation  of  the  actual  diffusion  length.  To  do  this,  the  beam 
source  was  modeled  as  a  Lorentz  function  while  the  luminescence  distribution 
was  modeled  as  an  exponential  function.  However,  the  model  was  a  first-attempt 
effort,  which  used  an  approximation  of  the  y  parameter  for  the  Lorentz  function 
representing  the  beam  source.  This  approximation  was  based  on  the  assumed 
value  of  45  cmVs'1  for  the  hole  carrier  mobility.  While  the  model  worked  and 
gave  a  better  fit  of  the  electron  carrier  mobility  values,  there  is  a  certain  degree  of 
inaccuracy  involved. 

To  address  the  problem,  there  is  a  need  to  determine  the  true  value  of  the 
y  parameter,  which  characterized  the  beam  source.  One  approach  is  to  use 
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sample  materials  with  known  transport  parameters  (i.e.,  mobility  and  lifetime)  for 
the  testing  of  the  model.  The  experimental  diffusion  length  is  a  convolved 
product  of  the  actual  diffusion  length  and  the  Lorentz  distribution.  Since  the 
actual  diffusion  length  is  known  (as  it  can  be  calculated  as  a  mobility-lifetime 
product),  different  values  for  the  y  parameter  can  be  fitted  to  see  which  gives 
the  closest  fit  to  the  experimental  result.  Once  the  value  of  the  y  parameter  is 
determined,  a  direct  deconvolution  of  the  coupled  signals  can  then  be  done 
instead  of  performing  the  two-stage  process  used  in  this  work. 
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